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We report on dependence of conductance and tunnelling magnetoresistance on bias voltage at 
different temperatures down to 2K in Co|Al2O3(10A)|Si(<5)|Al2O3(2A)|Permalloy magnetic tunnel 
junctions. Complementary low frequency noise measurements are used to understand the conduc- 
tance results. The obtained data indicate the breakdown of the Coulomb blockade for thickness of 
the asymmetric silicon layer exceeding 1.2A. The crossover in the conductance, the dependence of 
the tunnelling magnetoresistance with the bias voltage and the noise below 80K correspond to 1 
monolayer coverage. Interestingly, the zero bias magnetoresistance remains nearly unaffected by the 
presence of the silicon layer. The proposed model uses Larkin-Matveev approximation of tunnelling 
through a single impurity layer generalized to 3D case and takes into account the variation of the 
barrier shape with the bias voltage. The main difference is the localization of all the impurity levels 
within a single atomic layer. In the high thickness case, up to 1.8A, we have introduced a phe- 
nomenological parameter, which reflects the number of single levels on the total density of silicon 
atoms. 



PACS numbers: 

I. INTRODUCTION 

The discovery of large tunnelling magnetoresistance 
(TMR) at room temperature P, Q has strongly renewed 
the interest in spin tunnelling phenomena. Up to very 
recently the main efforts were concentrated on the in- 
crease of tunnelling magnetoresistance values by using 
ferromagnetic electrodes with the highest possible spin 
polarization (half metallic ferromagnetism) , searching for 
new types of insulating barriers (including the so called 
spin filters [H), or a combination of both approaches, 
where the ferromagnetic/insulator interface design could 
also play an important role for spin polarized tunnelling. 
The last approach has recently provided an enormous 
progress in the TMR values at room temperatures. It has 
been demonstrated that in epitaxial Fe/MgO/Fe mag- 
netic tunnel junctions (MTJ), where there exist condi- 
tions for coherent propagation of specific spin orbitals 
from one ferromagnetic electrode to the other one, TMR 
has reached experimentally values up to 410% [4]. This 
fact, supported by theoretical calculations predicting a 
TMR of more than 1000% have further increased both 
the fundamental and technological interest in spin polar- 
ized tunnelling. 

Another possible research direction, which remains 
however poorly explored, is related to tunnelling in com- 
plex (hybrid) junctions. Indeed, the manipulation of the 
barrier by doping with magnetic or nonmagnetic impu- 
rities, or inserting magnetic, for example see reference 



Q, nonmagnetic @ (even superconducting [HQ) quan- 
tum dots (QD), would add a new degree of freedom to 
spin polarized tunnelling and strongly enhance the ver- 
satility of spintronic devices based on spin polarized tun- 
nelling. Tunnelling in such hybrid, but non magnetic 
devices, has been intensively studied during the last two 
decades and especially in single electron transistors where 
the gate electrode is attached to quantum dots separated 
by two barriers each with a metallic contact (emitter and 
collector) Q . 

Recent theoretical studies of hybrid spintronic devices 
with two ferromagnetic leads contacting single or dou- 
ble quantum dots have revealed plenty of new interest- 
ing phenomena related with the interplay between mag- 
netic tunnelling processes and spin/charge accumulation 
on quantum dots. From the experimental side, some 
few groups have demonstrated that Coulomb interaction 
may indeed play an important role not only in ferromag- 
netic granular systems Q (as occurs in the correspond- 
ing nonmagnetic analogs), but also in ferromagnetic sin- 
gle electron tunnelling devices constructed either from a 
2D electron gas [l(| or when a single metallic nanoparti- 
cle is contacted by ferromagnetic electrodes [llj. It has 
also been reported that in ultrasmall nanoparticles, with 
a diameter between 5-10nm, their electronic structures 
(i.e. quantum effects) may also influence spin polarized 
tunnelling [12]. Other interesting examples of spin polar- 
ized tunnelling in hybrid structures include ferromagnetic 
leads contacting carbon nanotubes [l3[ or a C60 molecule 
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Actually, from the technological point of view, it is eas- 
ier to attach ferromagnetic leads via tunnelling barriers 
to an array instead of to a single quantum dot. This 
method, despite the evident drawback due to some dis- 
tribution in QD sizes and the corresponding charging en- 
ergies, adds evident versatility to the design of the exper- 
iment, allowing continuous tuning between two different 
regimes: (i) weak doping regime where QDs are substi- 
tuted by impurities and (ii) strong doping regime. In 
the second regime one would expect the QD charging en- 
ergy to reduce continuously with the average dot size, 
allowing sequential electron transport through an other- 
wise blocked channel (Coulomb blockade) in addition to 
direct tunnelling. 

Therefore, the knowledge of different mechanisms af- 
fecting the conductance, and especially its dependence 
on the applied bias, could be an important instrument 
for a comparative analysis of noise and transport. We 
should, however, note that the bias dependence of con- 
ductivity and TMR in real magnetic tunnel junctions, is 
still poorly understood. The most accepted model (la ], 
which does not take into account the possible defects in- 
side the barrier, predicts an attenuation of the polariza- 
tion due to the decrease of the difference between the 
height of the barrier and the bias voltage. However, the 
observed reduction in conductivity usually exceeds the 
predicted effect. 

Two theories have been developed in order to explain 
the anomalous reduction in TMR. The first one explores 
the existence of impurity states, which reduce the spin 
current polarization and influence the conductance of the 
junctions at low bias [16] . The second one involves inelas- 
tic tunnelling as the main origin of unpolarized current 
[l7j . This point of view has been also supported by the 
measurement of inelastic tunnelling in magnetic tunnel 
junctions [l8[- Later on Ding et al [l^| have also de- 
tected a tunnelling magnetoresistance in vacuum based 
magnetic tunnel junctions with a reduced dependence on 
the bias voltage. These results indicate a possible in- 
fluence of the impurity states on the polarization of the 
tunnel current. 

For the nonmagnetic tunnel junctions with (nonmag- 
netic) nanoparticles inside the barrier the presence of a 
zero bias anomaly was first reported and explained by 
Giaever (20| . The main mechanism responsible for the 
appearance of the threshold voltage is Coulomb block- 
ade, which controls two steps tunnelling. In this model 
the threshold voltage is distributed from V to a max- 
imum voltage V s . In the doped tunnel junctions V s is 
given by the size of the doping particles, which deter- 
mines a charging energy, given by the capacitance of the 
particles. This provides some distribution in the popula- 
tion of electrons inside the particles. 

The two steps tunnelling in a magnetic tunnel junction 
has been later treated using the other method developed 
by Glazman and Matveev [2ll . |%|. This theory states 
that the current is defined, in each particle, by the tun- 



nelling rates from one of the electrodes to the central 
particle, and from the island to the other electrode. As 
soon as conductance and spin current polarization are 
modified, a modification of the dependence on voltage of 
the tunnelling magnetoresistance may be expected. In- 
deed, when the bias voltage is increased, the number of 
allowed two steps processes is also increased. 

Previously to this work, Jansen et al. [23[ studied 
MTJs with Si nanoparticles up to 1.8A introduced in 
a symmetric position. They observed gradual suppres- 
sion of tunnelling magnetoresistance. In the present work 
Si particles were introduced asymmetrically inside the 
barrier. While symmetric doping effectively separates 
the barrier into two parts with similar tunneling rates, 
the asymmetric doping is expected to affect weakly the 
largest tunneling rate minimizing the influence of the 
nonmagnetic Si doping (at least for relatively weak dop- 
ing levels). 

This work presents an experimen- 
tal study of electron transport in 
C3o(100A)/AlaO 3 (10A)/Si(5)/Al 2 Os(2A)/Py(100A) 
hybrid magnetic tunnel junctions, where the largest 
barrier is 5 times larger than the short one. We observed 
a continuous transition between the weak (Si impurities) 
and strong (array of Si quantum dots) regimes. In order 
to discriminate the different conductance regimes, we 
study the temperature and bias dependence of both the 
conductivity and TMR as a function of Si doping. 



II. EXPERIMENTAL DETAILS. 

Details of sample preparation have been published pre- 
viously [23] . For silicon doped samples, the tunnel barri- 
ers were deposited in two steps. After deposition of the 
underlying Co electrode, a first tunnel barrier was formed 
by deposition and subsequent oxidation of 10 A of Al. 
Subsequently, sub-monolayer amounts of Si (8 in the fol- 
lowing) were deposited on the AI2O3 surface, followed by 
a second Al layer deposition (2 A) and oxidation, result- 
ing in a "(5-doped" Al2C>3|Si|Al2 O3 tunnel barrier. After 
the deposition of the barrier 100 A of permalloy were de- 
posited, in order to form the second magnetic layer. All 
the samples thicknesses were measured using a quartz 
monitor. The quality of the Al layers and the AI2O3 has 
been previously tested in Ref . (24 . |25| . In the following 
the positive bias voltage corresponds to the application 
of the voltage in the top electrode, while the negative to 
application of the voltage in the bottom one. 

Measurements were performed using a computer con- 
trolled system [26], which allows to detect the dynamic 
resistance, the DC value of the current and the voltage, 
and the noise in the device under study. Biasing of the 
samples was done at a constant current, applied using a 
calibrated source. It also allowed to modulate the ap- 
plied current. A square waveform was used in order to 
detect the transfer function of the line and the dynamic 
resistance of the junctions. The voltage response of the 
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FIG. 1: Barrier parameters evaluated using the BDR model. 
The graph (a) shows the barrier width and height in the left 
and the right axis, respectively, as function of the Si thickness. 
In the graph (b) is shown the asymmetry also as function of 
the Si thickness. In the inset we present the typical depen- 
dence on voltage of the conductance at room temperature. 
The dashed line is a parabolic fit used to extract the param- 
eters plotted in the graph (a) and (b). 



devices was amplified by using DC coupled low noise am- 
plifiers. The amplified signal was recorded in an analog- 
digital converter. 

The measurement of the noise uses the same biasing 
technique and the same low noise amplifiers, which were 
placed in the top part of a cryostat. The pre-amplihcd 
signals are further amplified by additional low-noise am- 
plifiers (Stanford Research SR560). A spectrum analyzer 
SR780 calculates the cross-correlation spectrum of the 
voltage noise, containing thermal, shot and 1// contribu- 
tions. The obtained dynamic resistance allows to convert 
the voltage noise into the current noise. Extrinsic noise, 
introduced by the amplifiers and the current source, was 
removed by using the data extracted from a careful cal- 
ibration performed on resistors at low temperatures. In 
the experiments we measured 9 different samples, 2 of 
each Si thickness, except for (5=0.6Awhere we only chara- 
terized 1. 



III. CONDUCTANCE, ZERO BIAS ANOMALY 
AND TUNNELLING MAGNETORESISTANCE 

At room temperature the dependence of the conduc- 
tance on the bias voltage fits well to a parabolic func- 
tion (see inset to Fig. QJ>). The parabolic depen- 
dence of the conductance, being due to direct tunnelling, 
was explained within the Brinkman-Dynes-Rowell (BDR) 
model [27j, which describes the tunnelling conductance 
as a function of the parameters of a trapezoidal barrier: 
the width, the average height and the difference between 
the sides of the barrier, such difference being known as 
asymmetry. 

The results of the fits are plotted in Fig. [TJ showing 
in general an increase of the barrier width with Si dop- 
ing. However, the increase of the width does not corre- 
spond with the deposited Si thickness. As to the height 
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FIG. 2: Resistance of the MTJ's at low temperature, graph 
(a) and room temperature, graph (a), as a function of the Si 
thickness (5(A)). The lines in the graphs are guides to the 
eyes. 



of the barrier, it clearly diminishes when the silicon is 
introduced inside the barrier. At present we have no 
clear explanation for this reduction. One of the possible 
reasons could be some decrease of the work function of 
the aluminium oxide. This explanation, however, con- 
tradicts the observed variation of the barrier asymmetry 
(Fig. [TJ, which decreases with the thickness of the Si 
layer. The disagreement between the deposited thickness 
and the obtained parameters using the BDR model could 
be attributed either to the simplification made by the 
model, which assums a trapezoidal barrier, a parabolic 
band structure and the WKB approximation, or to de- 
fects inside the barrier that should diminish the barrier 
height. However the obtained results present evidence 
that the Si layer affects the barrier properties beyond 
the errors committed in the estimation procedure. 

The tunnelling resistance of the studied Si doped MTJs 
measured at T =2 K and T =300 K, in the zero bias limit, 
is plotted in Fig. [2l The observed enhancement of the 
resistance when the temperature is lowered rules out the 
presence of pinholes even for the highest silicon thick- 
ness dil H^l . There is an increase of the resistance for 
8 >1.5 A, both at room and low temperatures, whereas at 
lower concentrations the resistance trend is to decrease. 
This fact indicates a change of the regime in the conduc- 
tance, which will be further confirmed in the dependence 
of the conductance on the bias voltage and the TMR. 



A. Dependence of the conductance on the voltage 
at low temperatures 

While at room temperature conductance is a parabolic 
function of the bias, at low temperatures (below 100 K), 
the dependence of conductance on the bias changes sub- 
stantially. On the one hand, all MTJs at low biases 
(V<30 mV) present a peak in the resistance. Usually 
such a peak is called zero bias anomaly (ZBA). On 
the other hand, high bias conductance regimes (above 
100 mV) show a strong variation with silicon thickness. 

To show more clearly the qualitative change of conduc- 
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tance regime with Si doping we present the bias depen- 
dence of the normalized conductance 



ZBA(%) = 100 



fi(100mV - fl(OmV) 
i?(100mV) 



(1) 



The dependence of ZBA vs. Si thickness, plotted in 
Fig. HI shows that the resistance peak, being weakly 
dependent on S for low Si thickness, strongly increases for 
the high doping region. The crossover region corresponds 
to the thickness of approximately 5—1.2 A. 

The experimental data presented above may be under- 
stood within the two-step model as follows. If the forma- 
tion of silicon islands starts for a Si thickness 6=1.2 A, 
then the enhanced resistance peak (ZBA) could be at- 
tributed to the appearance of a new energy scale in the 
electron transport through the barrier, related with the 
finite electron capacitance of the silicon islands, being 
practically absent for the small Si thickness range. This 
hypothesis is supported by the plot of the normalized bias 
voltage dependence of the conductance, Fig. [3l Clearly 
in the conductance of the junctions with higher Si dop- 
ing levels there is a crossover at a certain voltage. At 
higher voltage the dependence of the conductance be- 
comes less pronounced, whereas the junctions with low 
doping maintain the same behavior, as expected in a non- 
doped tunnel junction with electron transport due to di- 
rect tunneling. 

The change of the conductance regime is also evident 
from comparison of the lower bias conductance for un- 
doped MTJs and those with highest Si doping (corre- 
sponding to (5=1.8 A). As can be seen in the (b) graph 
in Fig. [31 both curves show structure in G(V) at low 
bias close to V=30 mV. The similarity of these weak 
anomalies both for the undoped and doped MTJs in- 
dicate their common origin, most probably related with 
electron conduction mechanisms through the aluminium 
oxide barrier. 



We have observed that generally, for all MTJs studied, 
the low bias conductance varies linearly with tempera- 
ture at low temperatures (T<20 K). Within the Coulomb 
blockade model this could be attributed to a variation 
in the thermally activated population of electrons inside 
the islands [3(|, which determines the slope of the con- 
ductance at zero bias. In brief, the conductance at zero 
bias and low temperatures could be expressed as 



G(0,T) <x / n(V ch )e 



~dV ( 



Ch 



n(V c 



Ch 



0)^T 

e 

(2) 

where eVch is the charging energy needed to intro- 
duce an electron in the metallic layer. At V =0 V the 
number of charged islands is given by the exponential 
term. Then, at low temperatures, the conductance is pro- 
portional to the number of charged particles in thermal 
equilibrium n (0). Although the temperature dependence 
of conductivity of the different MTJs studied varies in 
nearly two orders of magnitude, the normalized (to 2 K) 
low temperature slope in the linear dependence of con- 
ductivity vs. temperature was found to be much weakly 
dependent on the silicon thickness, with an average value 
of (5 ± 3) x 10~ 4 K _1 but with a rather large dispersion 
(Fig.©. 



B. Tunneling magnetoresistance 

The dependence of the tunnelling magnetoresistance 
on Si thickness is shown in figure [6l This plot represents 
the zero bias tunnelling magnetoresistance, obtained by 
using the following definition of TMR 
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FIG. 3: Dependence of the conductance on the bias voltage 
at T=2 K. The curves are normalized by the conductance 
measured at OmV in each curve. The plot (a) presents data 
up to high bias, where a mechanism related to the formation 
of Si islands becomes important (see the text). The graph (b) 
presents the low bias behavior, showing the presence of the 
same zero bias anomaly at V<30mV. 
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FIG. 4: Dependence of the quantity ZBA (defined in the 
text) on silicon thickness at T=2 K in the parallel (P) and 
the antiparallel (AP) magnetic state. The increased ZBA at 
higher doping layers reflects a different behavior, related with 
the formation of Si islands. 
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FIG. 5: The graph (a) plots the slope of the temperature 
dependence of the conductance normalized by its value at the 
lowest temperature. The graph (b) shows the low temperature 
linear dependence of the conductance on temperature. 
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FIG. 6: Dependence of the tunnelling magnetoresistance on 
the silicon thickness (5(A)) at three different temperatures. 
While presence of silicon reduces significantly the polarization 
at room temperature, decreasing the value of the TMR ratio, 
for temperatures below 80K the presence of the silicon layer 
seems not to affect the spin polarization, for 5 <1.2A: Some 
small influence appears at 8 > 1.2 A. 



TMR(%) = 100 



Rap — Rp 
R~p 



(3) 



where Rap and Rp are the resistance in the antiparallel 
and the parallel states, respectively. 

We have analyzed zero bias TMR vs. Si thickness for 
three different temperatures (300 K, 80 K, 2 K). The 
influence of the silicon doping on TMR is strongest at 
room temperature, suppressing TMR in nearly one order 
of magnitude for the highest Si thickness (1.8 A). The low 
temperature TMR values, however, remain nearly unaf- 
fected by the silicon. The step-like reduction of TMR 
at low temperatures (of about 10% at 77 K and 2 K) 
was observed for a Si thickness of 8 >1.2 A. This appar- 
ent reduction of the tunnelling magnetoresistance may be 
directly linked to variation of the ZBA with Si content 
(shown in Fig. 



Indeed, for small Si thickness, only direct tunnelling is 
possible, due to Coulomb blockade. This weakly changes 
TMR at low temperatures for small bias, which does not 
activate possible spin mixing due to the spin flip processes 
introduced by the Si. The Coulomb blockade is sup- 
pressed for the doping range of 8 >1.5 A, as indicated by 
conductance vs. voltage measurements (Fig. [3]), opening 
new conductance channels related to two-step tunnelling 
via the array of Si dots. The newly opened conductance 
channels create also a source of unpolarized current due 
to spin mixing and loose of spin memory of the electrons 
tunnelling through the array of Si dots. Suppression of 
Coulomb blockade just for the Si thickness bigger than 
1.2 A could be due to activation of a segregation process 
of Si atoms to nanometer scale dots, acting as a real ca- 
pacitance. This is contrary to the behaviour in low Si 
thickness junctions, where seems reasonable to suppose 
that Si could be more homogenously diluted inside the 
AI2O3 barrier in form of impurities and defects. In fact, 
two step tunnelling could, in principle, affect the con- 
ductance for both regimes discussed above. This is due 
to the unavoidable presence of defects inside the barrier, 
even without Si doping [31], which mix the spin currents 
and are a source of an unpolarized current. Therefore, 
this implies the presence of a finite characteristic spin-flip 
time on the defects and on the silicon layer both for the 
low and the high Si thickness regimes |32j. 

The analysis of bias dependence of TMR which follows 
further supports our hypothesis. In order to analyze the 
bias dependence of TMR as a function of Si thickness, 
we have found the voltage needed to suppress a TMR 
to its half (zero bias) value (i.e. to TMR(0 V)/2). This 
parameter, called Vtmr/2 is shown in FigJTl plotted as 
a function of the silicon thickness. Evidently, there is 
a crossover from a nearly constant Vtmr/2 regime be- 
low 8 =1.2 A, to a strongly decreasing one Vtmr/2 as 
a function of S, for Si thicknesses above one monolayer. 
The low Si doping regime with nearly constant Vtmr/2 
proves the presence of two step tunnelling through local- 
ized states, with character and density of states nearly 
unchanged up to Si thickness of 1.2 A. 

When 8 > 1.2 seems that the effective capacitance 
corresponding to localized states inside the barrier is re- 
duced, increasing dramatically the number of states avail- 
able inside the barrier for tunnelling. The new transport 
channels serve as a source of unpolarized current, which 
explains the much stronger voltage dependence of TMR 
for large Si doping. 



IV. THEORETICAL MODEL 

To account for the above discussed experimental fea- 
tures of electronic transport in tunnel junctions, we calcu- 
late now theoretically transport characteristics as a func- 
tion of bias voltage V using the model of rectangular tun- 
neling barrier with a thin layer of impurities (Si atoms) 
which create a number of impurity levels inside the bar- 
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FIG. 7: Voltage needed to reduce the tunnelling magnetore- 
sistance to 1/2 of its value at zero bias at T~2K. This value 
quantifies the dependence of the magnetoresistance on the 
bias voltage. A low value of Vtmr/2 means a strong depen- 
dence on bias voltage, while a high value means a weak one. 




We observe a diminished value of V 
starts to increase (see Fig. [4]). 



TMR/2 



when the ZBA 



rier. In the limit of low density of Si levels, the current 
is due to direct tunneling in each spin channel, and the 
current density can be calculated from the formula 



A. Impurity-mediated tunneling 

As the density of impurities grows, an additional mech- 
anism of tunneling through the impurity levels inside the 
barrier becomes more effective than the direct tunnel- 
ing. In the frame of the Larkin-Matveev [22] model, the 
resonant tunneling through the structure is described by 
the transition probabilities w^p. More specifically, w kp 
is the probability of transition from the state described 
by the wavevector p on the left side of the barrier to the 
state corresponding to the wavevector k on the right side 
(from now on we drop the spin index a referring to the 
tunneling in different spin channels) , and is given by the 
formula 



2tt 

Wkp = - 



£ 



£i + iTi 



S{e p - £ k ), (5) 



where T p i and T lk are the matrix elements for transitions 
between the states of the corresponding leads and of the 
i-th impurity, whereas F; is the width of the impurity 
level associated with tunneling from the localized level 
through the barrier. The sum in Eq. ([5]) runs over all 
impurities. The i factor is the imaginary unit. 

An important point is that the matrix elements T pl 
and Tit include a phase factor depending on the location 
of the impurity inside the barrier, 



i 



ki 



S l/2 



,-ikrRi e k z (zi-L/2) 



(2ttk) 



1/2 



(6) 



£ / de 



[2ra( £ -V,„)] 1 / ! /fi 



4:TT 2 h , /,, 

a u 

£ - h 2 kf/2m - V rc 
e - h 2 k 2 /2m - Vic 



1/2 
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where k; is the in-plane wavevector component of an 
electron incident on the barrier, t ka is the transmis- 
sion amplitude for an electron with wavevector k and 
spin a, while Vi a and V ra are the spin-dependent en- 
ergy band edges on the left and right sides of the junc- 
tion, respectively, which depend on the applied voltage 
as Vi a = V la0 + eV/2 and V r = V rM - eV/2 (here V /ct0 
and V ra o are the corresponding band edges at zero volt- 
age). Apart from this, the integration in Eq. ([4]) is over 
the electron energy e, and f(e) is the Fermi-Dirac distri- 
bution function. 

To calculate the tunneling probability for a nonzero 
voltage applied to the system we use a semiclassical ap- 
proximation for the wavefunction inside the barrier with 
a slope of potential. This is justified in case when the 
variation of the barrier height Vq is small at the electron 
wavelength A ~ h/(2mVo) 1 ^ 2 , which restricts the bias 
voltage to \eV\ <C m}l 2 V^ 2 L /h, where L is the barrier 
width. 



S l/2 



e i Pi .R, e - Pz (z l+ L/2) ^(2^)1/2^ 



(7) 



where S is the junction area, Ri and z\ are the in-plane 
and out-of-plane components of the i-th impurity posi- 
tion, k = (k/,ifc z ) (and similarly for p), while n is the 
inverse localization length of the impurity wavefunction. 
We assume that the impurities are randomly distributed 
in the plane, i.e. R; is a random variable, whereas z\ is 
the same for all impurities. 

Assuming that the energy level £; and the level width 
Tj do not depend on the position Rj, and averaging over 
Ri in the plane, we obtain the following formula for the 
electric current: 



i(v) =— yy^ iTl " Ti 

<7 k,p 



m) 2 + I? 



x£(£k-£p)[/(£k)-/(£k-ey)] 



(8) 



After calculating the average TkiXlpI , one finds that 
the current density j in the resonance-impurity channel 
consists of two terms, and can be written as 



, (£ k -£i) 2 + rf 

k,p 1 

xc5(£ k - £ P ) [/(e k ) 



Pi 



/(ek - eV)} 



(9) 
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FIG. 8: Total conductance per unit area, GjS, calculated 
as a function of the bias voltage for indicated areal impurity 
concentrations. The other parameters are described in the 
text. 

The first term in Eq. © is linear in the 2D impurity den- 
sity n and describes the transitions through completely 
isolated single levels. Such transitions do not conserve 
the in-plane components of p and k. The second term 
in Eq.© is nonlinear in n and describes the electron 
transitions through the impurity plane. For such transi- 
tions the corresponding in-plane components of electron 
wavevectors are conserved. 

In our calculations we include all three channels. The 
total conductance per unit area, G/S, is presented in 
Fig. [5] for parallel magnetic configuration and for indi- 
cated impurity concentrations. The width of the tun- 
nel barrier is taken as L = 1.2 nm, and the Si atoms 
are located within the plane of z\ = —0.4 nm, mea- 
sured from the center of the barrier. The energy struc- 
ture corresponds to the majority and minority bands in 
Co, Epi — 4.5 eV and Ep\ = 0.66 eV, respectively. 
The height of the barrier is assumed to be Vq — 1 eV. 
In turn, in Fig. [5] we show the reduced conductance, 
G(V)/G{V = 0). 

It should be noted that the conductance as well as 
the reduced conductance are slightly asymmetric with re- 
spect to bias reversal, see Fig. [5] and Fig. [5] This asym- 
metry is here associated with tunneling through single 
impurity levels which are located asymmetrically within 
the barrier. In our calculations we assumed that the im- 
purities are located in the barrier close to the interface 
between the barrier and one of the electrodes - like in ex- 
periments discussed above. Such an asymmetry of impu- 
rity position with respect to the center of tunnel barrier 
leads to the asymmetry of the conductance G(V). 



B. Role of Coulomb interaction 

When the density of impurity levels n grows, the con- 
ductance calculated within the model described above 
and including both direct and impurity-mediated tunnel- 
ing increases monotonically (see Fig. [SJ . The experimen- 
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FIG. 9: Reduced conductance, G(V)/G(0), as a function of 
bias voltage V for indicated density of Si atoms. The other 
parameters as in Fig. [S] 



tal data (Fig. [2]), however, reveal a rather sharp increase 
of the resistance to much higher values when the den- 
sity n crosses a critical value n cr corresponding to nearly 
complete filling of one atomic layer with Si atoms. The 
results for TMR presented in Figs. [6] and [7] also show that 
the physics of tunneling is substantially different for the 
density of Si atoms corresponding to complete filling of 
the plane. 

We assume that the physical reason of such transition 
is related to a dramatic increase of the role of Coulomb 
interaction in the conductance through the Si levels. This 
effect can be also described by the decrease of the num- 
ber of active levels, which are able to transmit electrons 
through the barrier. Indeed, if some Si atoms form a 
small cluster, then the cluster of several atoms acts as 
a single level for transmission because the Coulomb in- 
teraction prevents two or more electrons to occupy the 
same cluster. Thus we can assume that Eq. © describes 
the conductance as a function of the density of effective 
levels, n — > n e /, corresponding to the number of im- 
purity clusters. As the density of Si atoms approaches 
the critical value n cr , the value of n e / decreases rapidly. 
If n > n cr and the density of Si atoms keeps growing 
then it corresponds to increasing thickness of the layer 
completely filled with the Si atoms. In such a case the 
Coulomb interaction is suppressed as there are no small 
clusters anymore, and the conductance can be described 
using a model of three-layer structure with well-defined 
properties of each of the layers. One can expect that the 
properties of the Si layer in such a structure are close to 
those of layered amorphous Si. 

It should be emphasized that the direct tunneling is 
also suppressed in the vicinity of n ~ n cr . This is re- 
lated to the Coulomb repulsion of electrons transmitted 
through the barrier from the charged impurity clusters, 
so that the electrons can tunnel through the barrier only 
in those areas, which are free from the impurity clusters 
or islands. One can describe this by a local increase of the 
tunneling barrier in the areas filled with clusters. This 
effectively leads to an increase of the average tunneling 




0.0 0.05 0.1 0.15 

6 (nm) 

FIG. 10: Variation of the resistance per unit square as a func- 
tion of Si layer thickness. 

barrier. 

Using the above described ideas we have calculated the 
resistance as a function of 8, taking into account all three 
channels of conductivity as described above, but instead 
of the number of impurities n we put in Eqs. ([H]),© an ef- 
fective number of levels n e f , which we assumed to change 
rapidly from n e j = n at n < n cr ~ 1.1 x 10 14 cm~ 2 
to a constant value 10 13 cm -2 . We also corrected the 
contribution due to direct tunneling making it strongly 
dependent on n in the vicinity of the transition point 
ri = n cr . More specifically, we reduced this contribution 
for n > n cr by modeling the dependence of the tunnel 
barrier Vq on n: for n < n cr we take Vo = const inde- 
pendent of n but for n > n cr we assume that this value 
increases by 0.6 eV, which corresponds to suppression of 
the direct tunneling due to the 'screening' from the large 
impurity clusters within the Si layer. The results for the 
resistance as a function of S are presented in Fig. [TUJ As 
one can note, the theoretical curve is qualitatively similar 
to the experimental one. 

Further improvement of the model can be made by 
taking into account the dependence of the impurity den- 
sity of states v on energy in the vicinity of the Fermi 
level, e = a. One can assume that the function v{e) 
has a minimum near e — n in accordance with the 
shape of the density of states (DOS) in amorphous Si 
(see, for example, Refs. [H, |34|). The DOS profile 
in the vicinity of the minimum at e = fx in amor- 
phous a-Si can be approximately presented as a sum 
of the DOS tails related to the conduction and valence 
bands, u(e) ~ v Q ( e -(e-^+A)/ e „ + e (e-M-A)/ £o ) ; w here 

A ~ 0.25 eV -1 atom _1 , and £o — 100 meV. In the vicin- 
ity of the minimum at e = one can use the parabolic 
approximation 

K^)^Km)(i + ^p) (io) 

where v(fi) = 2isoe~ A / eo . 

After substituting the constant DOS v{p) by the ap- 
proximate function (JTUJ) with Eq — 100 meV, we finally 
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FIG. 11: Reduced conductance in the model with Coulomb 
repulsion and the DOS of amorphous Si layer taken into ac- 
count. 

find the dependence presented in Fig.QT] The calculated 
dependence shown in Fig.QT] is m reasonable qualitative 
agreement with the experimental curves Fig. [31a. 



V. DISCUSSION 

Several attempts have been made previously trying to 
understand the possible role of Coulomb blockade in the 
tunnelling current and magnetoresistance for ferromag- 
netic leads contacting a quantum dot [H, Ha|- These 
works predict an oscillation in TMR(V) with a period 
given by the charging voltage. Experimentally, for tun- 
nelling through an array of dots, the oscillatory behaviour 
has been reported only for the magnetic tunnelling junc- 
tions with a barrier doped with cobalt nanoparticles. 
Those junctions had rather small ferromagnetic elec- 
trodes area A=0.5x0.5 /im 2 with the conduction almost 
completely blocked at low bias [9|, and TMR oscillat- 
ing with the period predicted theoretically. Other ex- 
periments, which also studied spin dependent electron 
transport through arrays of dots doping the barrier of 
magnetic tunnel junctions, have reached conclusions in 
respect to conductance similar to those reported here, i.e 
observation of a non-oscillating increase of the conduc- 
tivity with bias, when Coulomb blockade is suppressed 
by the applied voltage [H, 0, [n| 

The above cited reports employed different devices 
in order to study spin dependent tunnelling through a 
medium controlled by Coulomb blockade, and two of 
them specifically Ref. [5j and Ref. [9( used qualitatively 
similar magnetic tunnel junction devices. In these papers 
a granular film, consisting of nanometer size cobalt par- 
ticles (with radius close to ^2.5 nm) was embedded in a 
matrix of aluminium oxide. This array of Co dots was 
deposited on top of the aluminium oxide barrier (2.7nm 
in the first case [f| and 1-2 nm in the second one Q)and 
was covered by the second aluminium oxide barrier. The 
top barrier was different in the studies mentioned above. 
While the first paper [5j used a ~1.5 nm thick secondary 
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FIG. 12: Typical random telegraph noise process. The plot 
(a) shows the power spectrum of the process. It is dominated 
by a Lorentzian added to a 1/f noise background (straight 
line in the graph). The graph (b) of the figure shows a typical 
time series. The two states fluctuation is clearly seen in the 
histogram (shown in (c)). The histogram shows the number 
of counts at a certain bin of voltage. Two peaks correspond 
to two states of different conductance. 



barrier, the later work did not use any alumina barrier to 
cover the Co nanoparticles, which probably produced an 
uncontrolled secondary barrier. Another difference is re- 
lated to the junctions area. In the first study a rather big 
(4.5 xlCP 2 mm 2 ) area MTJ was used, while the second 
work studied Coulomb blockade controlled spin depen- 
dent tunnelling in MTJ's with much smaller area junc- 
tion. All different studies of the spin dependent trans- 
port in MTJ's with nanoparticle doped barrier, with the 
exception of Yakushiji et al., reported a staircase de- 
pendence of the IVs presumably due to single electron 
charging effects, and did not show any conductance os- 
cillation. 

Although our samples have a silicon S layer inside the 
barrier, instead of a magnetic S layer introduced in the 
previous reports, the observed behaviour of the conduc- 
tance vs. voltage is in general similar to the data reported 
for spin tunnelling through arrays of magnetic nanoparti- 
cles Q , with exception of the dependence of TMR on bias 
voltage which is smoother (Vtmr/2 ^0.5 V, [f|) than in 
our samples. 

A further confirmation of the role of the Coulomb 
blockade in our samples could be observed in the noise at 
low frequency in the studied samples. Whereas for low 
Si doping (6 < 1.2 A), the power spectrum at low tem- 
peratures is "white" (i.e. nearly frequency independent) 
and corresponds well to the shot noise expected for direct 
tunnelling in a tunnel junction (or for two-step tunnelling 
through strongly asymmetric barriers), for S >1.5 Aa 
random telegraph noise (RTN), Fig. contribution be- 
comes evident for bias voltages above a critical value. 
The appearance of RTN for 5 >1.2 Amight be under- 
stood as a consequence of the suppression of Coulomb 
blockade. 



RTN has been previously reported for nonmagnetic 
tunnel junctions and some other devices such as field ef- 
fect transistors or quantum dots connected to metallic 
leads [38]. As to the tunnel junctions, RTN has been 
usually attributed to resistance fluctuations due to a sin- 
gle or few fluctuators (3^. The RTN was usually found 
for rather small area (below 1/im 2 ) junctions and low 
temperatures, because in this case the tunnel resistance 
is controlled by a few fluctuating defects, providing two- 
state fluctuations of the resistance. In the case of rather 
large tunnel junctions, as in the present study, the ob- 
servation of RTN could not be described by the above 
models, involving direct influence of single or few defects 
fluctuations on the resistance. A non-uniform current 
distribution, induced by pin-holes, which are a source of 
"hot spots" just before the MTJs are broken down by the 
intensity of current, would neither explain the observed 
voltage dependence of the RTN. Indeed, our experimental 
data, particularly the current-voltage characteristics and 
the temperature dependence of the conductance show ab- 
sence of pinholes and the above mentioned "hot spots". 

As discussed above, for thick enough silicon layers 5 > 
1.2 A the effective capacitance of the Si dots becomes 
small enough to break down the Coulomb blockade above 
a certain bias voltage. This increases the electron popu- 
lation of the island due to the two steps tunnelling events 
and enhances the tunnelling conductance. 

It is evident that in the system under study the ca- 
pacitance of the Si dots and, correspondingly, the Si dots 
population, should be distributed over the MTJ area, 
providing a possible variation of the local tunnelling cur- 
rent as a function of the spatial coordinate. In addition, 
the two-level systems situated close to the Si dots seem 
to introduce also time dependent fluctuations or RTN 
in the tunnelling current through these dots. The un- 
avoidable dependence of the tunnelling current on the 
coordinate may enhance the effective contribution to the 
overall conductance from only a few fluctuators, resulting 
in a noise contribution additional to 1/f due to effective 
"amplification" of some conductance fluctuations respon- 
sible for RTN from these few defects. This is represented 
schematically in the Fig. Q2] the change in the charge 
in the defect yields two states with different tunneling 
rates (r), hence different conductance. The origin of the 
different times lies in the different levels inside the dots, 
which are determined by the capacitance, thus by the lo- 
cal environment. Of course, in our tunnel junction there 
are many defects, but at low temperatures there are only 
several of them active because the trapping-detrapping 
process uses to be thermally activated. 



VI. CONCLUSIONS 

To conclude, we have carried out an extensive study of 
electron transport in Co| AI2O3IPV magnetic tunnel junc- 
tions asymmetrically doped with Si. Our experimental 
data suggest that the observed behaviours of conduc- 
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electron population, even though their capacitances are 
characterized by rather narrow size distributions, giving 
rise to some distribution of the maximum threshold volt- 
age suppressing Coulomb blockade ±V c h- This leads to 
the observed zero bias anomaly in our samples. On the 
other hand, the variation of the electron population in 
each particle may explain the suppression of the pre- 
dicted [35] Coulomb oscillations, which could be present 
only for a constant equilibrium electron population in Si 
particles across the junction area. One of the ways to 
reach a more uniform electron population in the Si parti- 
cles could be the reduction of the area of the junctions as 
in Ref. [9j, implying tunnelling through a smaller array 
of nanoparticles with a size distribution narrower than in 
the present case. 



FIG. 13: Schematic explanation of the noise observed at low 
temperatures. The presence of traps in our barrier modifies 
the levels inside our dot. This effect leads to two resistance 
states depending on what is the state of the defect. This effect 
in visible at low temperatures because only few processes are 
activated, hence the noise is RTN type as is clearly seen in 

Fig. m 

tance vs Si doping are closely related, clearly indicat- 
ing a suppression of the Coulomb blockade regime for Si 
layer thicknesses above about one monolayer. Although 
no staircase behavior of the conductance was observed 
for the regime of suppressed Coulomb blockade, as in 
some other systems [9|, the observed behavior of tun- 
nelling conductance is rather similar to the one reported 
long time ago by Giaever [20l | , which was successfully ex- 
plained as due to the presence of a large amount of par- 
ticles inside the barrier. For a fixed bias voltage applied 
in equilibrium, these particles might have very different 
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